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ABSTRACT 


The  Infrared  absorption  of  water  vapor  has  been  calculated 
over  a  wide  r  nge  of  path  length,  pressure,  and  temperature  from 
1000  to  10, f  .0  cm**^.  The  results  are  presented  In  extensive  tables 
which  give  the  transmission  over  Intervals  of  2.5  cm~^  at  wave  numbers 
below  3400  cm"^.  The  average  transmission  over  Intervals  of  20,  50, 
and  100  cm'^  Is  given  over  the  complete  frequency  range.  The  cal¬ 
culated  data  representing  the  Intensity,  position,  and  number  of 
spectral  lines  In  a  given  frequency  Interval  Is  tabulated  In  a  form 
v^lch  Is  useful  for  transmission  calculations.  All  lines  from  all 
Isotopic  species  having  an  Intensity  greater  than  10~^  of  the  strongest 
line  In  an  absorption  region  were  Included  In  the  calculation.  The 
quasi-random  model  for  the  representation  of  band  spectra  Is  discussed 
In  detail.  Among  Its  features  are  the  possibility  of  representing 
the  absorption  by  a  series  of  lines  vdilch  are  arranged  neither  at 
random  nor  In  a  regular  way. 
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SECTION  1 


INTRODUCTION 


The  amount  of  Infrared  radiation  which  ia  transmitted 
through  the  earth's  atmosphere  is  determined  largely  by  the  be¬ 
havior  of  two  molecules:  water  vapor  and  carbon  dioxide.  In  this 
report  a  method  by  which  the  spectral  transmission  of  water  vapor 
has  been  calculated  over  a  wide  range  of  Infrared  frequencies  is 
described.  The  results  which  were  calculated  on  high  speed  digital 
computers  are  presented  in  the  form  of  extensive  tables.  A  coaq>anion 
report  (SSD-TDR-62-127  -  Volume  III)  presents  similar  data  for  carbon 
dioxide. 


The  most  extensive  transmission  tables  which  have  been 
previously  calculated  are  those  of  Ballard,  Larmore,  and  Passman.^ 
Their  tables  cover  a  wavelength  range  from  0.3  to  7 .0  fJL  at  intervals 
of  O.lju  .  The  results  are  for  one  temperature  (300°K)  and  one 
pressure  (1  atmosphere).  They  assumed  that  the  transmission  is  given 
by  the  Elsasser  model.  In  particular,  they  used  the  error  function 
approximation  for  the  Elsasser  model  which  is  of  limited  accuracy 
at  the  higher  pressures  and  shorter  path  lengths.  Furthermore,  the 
Elsasser  model  itself  is  a  poor  representation  of  the  water  vapor 
spectrum  because  the  individual  lines  are  not  arranged  with  the 
regularity  required  by  this  model. 

Tables  of  water  vapor  transmissions  averaged  over  20  cm*^, 
SO  cm~^,  and  100  cm~^  are  presented  respectively  in  three  groups. 
Entries  for  the  20  cm~^  averages  are  presented  every  10  cm~^  from 
1000  cm~^  to  10,000  cm'l.  nie  50  cm~^  and  100  cm~^  average  entries 


^S.  S.  Ballard,  L.  Larmore,  S.  Passman,  Fundamentals  of  Infrared  for 
Military  Application,  The  Rand  Corporation,  R-297,  March  31,  1956  (C). 
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are  presented  every  25  cm"l  and  50  cm“l  respectively  in  the  same 
frequency  range.  These  tables  contain  results  for  three  different 
temperatures  (300°K,  2500K,  and  200OK)  with  entries  for  each 
temperature  at  seven  pressures  and  fifteen  amounts  of  H2O.  The 
range  of  temperatures,  pressures  and  path  lengths  cover  the  range 
of  values  which  occur  in  the  earth's  atmosphere. 

Tables  of  the  transmission  at  intervals  of  2.5  cm~^  from 
1000  cm'l-  to  3400  cm'l  are  also  presented.  These  calculations  were 
carried  out  using  the  quasi-random  model  which  is  capable  of  repre¬ 
senting  the  spectral  transmission  of  molecular  gases  quite  accurately. 
A  detailed  description  of  this  model  is  given  in  Section  3. 

In  Section  2  the  molecular  structure  of  H2O  is  described 
together  with  the  rotational  and  vibrational  motions.  The  intensi¬ 
ties  and  frequencies  of  the  spectral  lines  are  described. 
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SECTION  2 


MOLECULAR  STRUCTURE  OF  WATER  VAPOR 


2.1  VIBRATIONAL  ANALYSIS 

The  water  molecule  Is  a  non-linear  molecule  with  three 
non-degenerate  elgen-vlbratlons,  two  of  which  are  symmetric  and  the 
other  antisymmetric.^  The  vibrational  potential  constants  have 
been  measured  to  second  order.  The  vibrational  energy  to  this  order 
is  given  by 


i  3 

L  C0.(v  4%)  +  T.  E  Xn(v  4%)(v  +%) 
i=l  ^  i-1  j  >i  ^  J  ^ 


where  the  constants  0)^^,  x^.  are  listed  in  Table  1.  is  the 
vibrational  energy  in  cm"l  and  v^  are  the  vibrational  quantum 
numbers.  The  data  for  the  main  isotopic  species  HHO^^  (abundance 
0.998)  were  obtained  by  Darling  and  Dennison. ^ 


^G.  Herzberg,  Infrared  and  Raman  Spectra  (D.  Von  Nostrand  Company, 
Inc.,  New  York,  1945). 

^B.  T.  Darling  and  D.  M.  Dennison,  Phys.  Rev.  128  (1940). 
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TABLE  1 


VIBRATIONAL  CONSTANTS 


Isotope 

HHO^^ 

HHO^® 

HHO^^ 

HDO^^ 

60 

3825.32 

3816.82 

3820.52 

2824.32 

0)2 

1653.91 

1647.21 

1650.41 

1440.21 

60  3 

3935.59 

3919.69 

3927. 19 

3889.84 

^11 

-  43.8 

-  43.6 

-  43.7 

-  43.4 

-  20.0 

-  19.9 

-  19.9 

-  8.6 

^13 

-155.0 

-154.0 

-154.5 

-  13.1 

X22 

-  19.5 

-  19.3 

-  19.4 

-  11.8 

X23 

-  19.8 

-  19.6 

-  19.7 

-  20.1 

*33 

-  46.8 

-  45.9 

-  46.1 

-  82.9 

18  17 

For  the  species  HHO  (abundance  .002)  and  HHO  (abundance 
the  constants  were  calculated  directly  from  the  relations^ 

.0004) 

1 

f  1 

2  sln^a) 

8.520  X  10^ 

(2) 

'^2  2 

4  7T  C 

^m(H) 

m(0)  1 

2  2 

00  +  00^  = 

1  / 

14.690x10^-6. 140xl05slna  .  21. 232x10^ cos^ a 

47T^c^  \ 

m(H) 

T 

m(0) 

(3) 
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1 


(4) 


4^4 
16  77  c 


2m(Hl  +  m(0^ 
[m(H)]^  m(0) 


33.824  X  10" 


2 

cos  a 


and 

(isotope)  (isotope) 
(isotope)  ^  ^ _ 

y  y  coj^cDj 


(5) 


where  m(H)  and  m(0)  are  the  masses  (in  grams)  of  the  hydrogen  and 
oxygen  atoms  respectively  and  2a  ~  105^^  is  the  angle  between  the 
H  atoms  at  the  position  of  the  0  atom.  The  numerical  constants  in 
Eqs.  (2  -  4)  and  the  value  of  a  are  obtained  from  the  measured 
values  of  for  The  constants  for  the  species  HDO  (abundance 

3  X  10"4)  were  obtained  by  Benedict,  Gailar  and  Plyler.^  The  other 
isotopic  species  of  water  are  not  sufficiently  abundant  to  be 
considered  here. 


For  many  molecules  accidental  degeneracies  occur  when  two 
or  more  vibrational  states  have  the  same  symmetry  type  and  nearly 
the  same  energy  (see  Chapter  11.5c  of  reference  1  or  SSD-TDR-62-127 
-  Volume  III  ).  This  situation  does  not  occur  in  H2O  for  any  of  the 
fundamental  vibrational  states  but  does  occur  for  the  states  (002) 
and  (200)  and  for  all  subsequent  pairs  of  states  (vj^  -  2,  V2,  V3  +  2) 
and  (v]^,  V2,  V3).  For  the  CO2  molecule,  vdiich  is  considered  in 
SSD-TDR'62-127  -  Volume  III,  this  resonance  occurs  between  the  funda¬ 
mental  state  (100)  and  the  state  (020).  Furthermore  for  CO2  the 
perturbation  associated  with  the  resonance  is  produced  by  the  third 
order  terms  in  the  potential  idiile  for  H2O  the  first  contribution 
comes  from  the  fourth  order  terms.  Thus  the  effect  is  less  significant 
for  H2O  than  for  OO2.  Darling  and  Dennison^  have  shown  that  the  off- 
diagonal  elements  in  the  vibrational  energy  matrix  for  H2O  are  given  by 


W 


Vi,V2,V3 


Vi-2,V2,V3+2 


^  Vv^(v^-1)(v3+1)(v3+2) 


(6) 


where  7  ■  74.4  cm"^’.  For  the  species  HDO,  20)^^  j*  2CO3,  but 


W.  S.  Benedict,  N.  Gailar  and  E.  K.  Plyler,  J.  Chem.  Phys.  24, 
1139  (1956). 
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I 


20)2  »  00^,  thus  the  resonance  occurs 

V2  -  2,  V3),  (vj,  V2,  V3).  Benedict, 
that  for  HDO 


between  the  states 
Caller  and  Plyler^ 


(Vj  +  1, 

have  shown 


V,  ,v  ,v,  _ 

^Vj+1,V2-2,V2  “  2  "V  (V^+1)V2(V2-1) 


(7) 


where  7  =  26.1  cm 

The  elements  of  the  vibrational  energy  matrix  can  be 
evaluated  from  Eqs.  (1)  and  (6)  or  (7).  The  energies  of  the  vibra¬ 
tional  level  are  the  eigenvalues  of  this  matrix.  A  calculation  of 
the  vibrational  energies  has  been  programned  for  the  IBM  709  computer. 
The  vibrational  energies  have  been  obtained  in  this  manner  for  all 
four  Isotopic  species  up  to  the  level  (013).  Levels  of  higher  energy 
than  this  were  not  considered  because  for  the  temperatures  of  Interest 
(  S  300°K)  only  the  levels  (000)  and  (010)  are  significantly  populated. 
Furthermore,  the  vibrational  transitions  between  (010)  and  all  levels 
higher  than  (013)  correspond  to  frequencies  greater  than  10,000  cm'^ 
which  Is  beyond  the  frequency  range  considered  In  the  current  study. 


2.2  ROTATIONAL  ANALYSIS 


In  regard  to  rotation  the  water  molecule  Is  an  asymmetric 
top  molecule  since  Its  three  principal  moments  of  Inertia  are  different. 
For  such  a  molecule  the  rotational  energies  cannot  be  expressed  as 
simple  functions  of  the  rotational  quantum  numbers. Unlike  the 
symmetric  rotator  there  are  no  rotational  degeneracies  In  the  energy 
levels  of  the  asymnetrlc  rotator.  In  addition  there  Is  no  "good" 
quantum  number  comparable  to  the  component  of  angular  momentum  about 
the  molecular  axis  of  the  symmetric  top  molecule.  Thus  the  spectrum 
of  rotational  energy  levels  which  Is  associated  with  each  vibrational 
level  cannot  be  described  by  an  analytic  expression. 


For  a  given  asymmetric  molecule  an  equivalent  symmetric  top 
molecule  Is  defined  as  one  for  which  the  two  smaller  moments  of  Inertia 
are  equal.  This  moment  is  taken  to  be  the  arithmetic  average  of  the 
two  smaller  moments  of  the  asymmetric  molecule.  In  order  to  describe 
the  rotational  energy  spectrum  of  an  asymmetric  molecule  it  has  been 
customary  to  use  the  equivalent  symmetric  top  molecule,  with  Its  good 
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quantum  numbers  and  analytic  energy  expressions,  as  an  approximation.  * 

It  Is  clear  that  this  Is  permissible  only  If  the  asymmetry  Is  slight. 

For  H2O,  which  Is  quite  asymmetric,  the  equivalent  symmetric  top  approxi¬ 
mation  cannot  In  general  be  used  to  obtain  the  correct  rotational 
structure.  However,  In  the  present  study,  the  exact  details  of  the 
rotational  spectrum  are  not  needed.  The  detailed  arrangement  of  the 
spectral  lines  within  a  finite  frequency  Interval  Is  not  needed.  If  the 
number  and  Intensity  of  the  lines  In  the  Interval  Is  known.  Thus  the 
symmetric  top  approximation  has  an  application  here  for  the  calculation 
of  positions  and  Intensities  of  the  weaker  spectral  lines.  However, 

It  Is  clear  that  there  will  be  a  lower  limit  on  the  size  of  the  fre¬ 
quency  Interval  which  may  be  used.  This  point  Is  discussed  In  more 
detail  In  section  2.3. 


In  the  present  calculation  of  water  vapor  transmission  the 
actual  asymmetric  lines  were  used  for  the  strongest  lines  In  the  lower 
portion  of  the  frequency  range  considered.  The  equivalent  symnetrlc 
top  approximation  was  used  only  for  the  weak  lines  In  this  portion  of 
the  spectrum  and  for  the  frequency  range  above  3400  cm~^  where  theo¬ 
retical  calculations  of  the  actual  water  vapor  spectrum  are  not 
available.  This  aspect  of  the  calculations  Is  discussed  In  detail 
In  section  2.3. 


It  can  be  shown  that  for  a  symmetric  top  molecule  the  rota¬ 
tional  energy  Is  given  by^ 


F  (J,K)  -  i  J(J+1)  +  (B^  -i  )K^  (8) 

2  ^1  ^"2 

where  Fy,  By^,Bv2  are  all  measured  In  cm  J  Is  the  total  angular 

momentum  quantum  number  vdtlch  takes  on  the  values  K,  K  +  1,  ...,  and 
K  Is  the  component  of  angular  momentum  about  the  symmetry  axis.  If  the 
symmetric  top  is  used  as  an  approximation  for  an  asymmetric  molecule, 

then  the  rotational  constant  Bu  Is 

^2 


+  B  ) 
"^3 


(9) 


4 

P.  C.  Cross,  R.  M.  Halner  and  G.  W.  King,  J.  Chem.  Phys.  27  (1943); 
Ibid.  12,  210  (1944). 
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vrtiere  .  are 
given  by^ 


the  rotational  constants  of  the  as3nnroetrlc  molecule 


T. 

nm  m 

m*l 


(10) 


h  6 

vrtiere  Bg  ■  - 2 -  >  “  “  1>2,3  and  are  the  principal  moments  of 

Inertia  corresponding  to  the  equilibrium  position  of  the  vibrating 
molecule.  The  constant  describes  the  effect  of  variation  of  the 
moments  of  Inertia  caused  by  the  vibration.  The  Integer  n  refers  to 
the  3  principal  axes  of  rotation  and  m  refers  to  the  3  modes  of 
vibration. 


The  numerical  values  of  the  rotational  constants  have  been 
determined  by  Benedict,  Gallar  and  Plyler^  and  from  the  Isotopic 
formulas 


g  (Isotope) 


m 


2/3 


(11) 


and 


^(Isotope)  ^ 
nm 


a  (  p  )■ 

nm  n' 


(12) 


where 


00  (Isotope) 
n 


This  data  Is  given  In  Table  2. 
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TABLE  2 


ROTATI(»)AL  CONSTANTS 


Species 

Constant 

(cm“l) 

HHO^^ 

HHO^® 

HHO^^ 

HD0‘‘ 

B 

27.379 

27.1910 

27.2790 

22.6906 

^1 

B 

22.0550 

11.9722 

12.0109 

7.8398 

^2 

“ll 

0.750 

0.745 

0.747 

0.253 

“21 

-2.941 

-2.905 

-2.922 

-1.798 

“31 

1.253 

1.238 

1.245 

1.087 

®12 

+ 

0.2199 

0.2185 

0.2191 

0.1544 

“22 

+ 

-0.0104 

-0.0103 

-0.0104 

-0.3800 

“32 

+ 

0.1113 

0.1099 

0.1106 

0.0503 

+ 


O-  ■  +  *_  )»  which  corresponds  to  using  Eq.  (9)  directly. 

2.3  INFRARED  SPECTRUM 

The  infrared  spectrum  of  water  vapor  consists  of  spectral 
lines  corresponding  to  transitions  bettwen  various  vibration- rotation 
energy  levels.  The  frequency  of  a  given  line  0} ,  is  proportional  to  the 
difference  in  the  vibrational  plus  rotational  energy  of  the  final 
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and  initial  states.  The  intensity  of  a  line  can  be  characterized  by 
the  integrated  absorption  per  unit  pressure  given  by^ 


3 

STT 

N 

r  E  _i 

VR 

1  M  iM 

3hc 

0)  TT-  exp 
QP 

"  kT 

K  • 

1  M  1 

1  -  e«P(-  ^ 


(13) 


where  N  is  the  number  of  molecules  per  unit  volume,  Q  is  the  partition 
function,  p  is  the  pressure,  M  is  the  dipole  matrix  element  corre¬ 
sponding  to  the  transition,  and  Eyg  is  the  vibrational  plus  rotational 
energy  of  the  Initial  state. 


The  selection  rules  which  govern  allowed  transitions  for  a 
triatomlc  asyninetric  molecule,  such  as  H2O,  are  given  by  Herzberg.^ 

For  the  molecular  species  HHO  there  are  Just  two  types  of  rotational 
transitions  allowed.  These  correspond  to  the  so-called  parallel 
(  AK  ■  0)  and  perpendicular  (AK  >  t  1)  transitions  of  a  synmetrlc 
top  molecule.  Bands  made  up  of  transitions  of  one  of  these  types 
approach  as  limits  the  parallel  or  perpendicular  bands  of  the  equiva¬ 
lent  symmetric  top  molecule.  They  also  have  the  same  general  structure. 
For  an  asymmetric  molecule  such  as  HDO,  a  given  vibrational  band  will 
contain  rotational  transitions  of  both  types.  The  corresponding  band 
of  the  equivalent  sjrmmetric  top  molecule  will  be  a  hybrid  band  con¬ 
taining  both  AK  -  0  and  AK  -  t  1  transitions.  (The  vibrational  bands 
of  HHD  molecules  are  not  hybrid  bands.)  The  familiar  rotational  selec¬ 
tion  rule  AJ  ■  0,  1*  1,  applies  to  all  transitions. 


In  addition  to  the  transition  frequency  u)  >  the  matrix 
element  M  cannot  be  expressed  explicitly  in  terms  of  "good"  quantum 
numbers  for  an  asymnetric  top  molecule.  Thus  the  expression  for  the 
integrated  absorption,  Eq.  (13)  is  quite  untractable.  On  the  other 
hand,  Eq.  (13)  becomes  relatively  simple  for  the  symmetric  top  molecule. 
In  fact  for  the  teiiq>eratures  of  Interest  here  the  interaction  of  vibra¬ 
tion  and  rotation  may  be  neglected  in  the  matrix  element  M.  Further¬ 
more  at  infrared  frequencies,  exp  [  .  _]<<!.  Taking  these  slmpllfl- 

kT 

cations  into  account,  Eq.  (13)  for  a  symmetric  top  molecule  reduces  to 


,v'K' J' 
*vKJ 


^v'K’J' 

vKJ 


V  V 

CO 

V 


exp 


fhcF  (J,K)] 

[  -v-j 


(14) 


v'  K'  J'  v' 

where  (U  1’,.  is  the  frequency  of  the  line,  O)  is  the  vibrational 

V  f  K|  w  V 


^S.  S.  Fenner  and  A.  Thomson,  Tech.  Rep.  No.  25,  Contract  No.  Nonr-220(03) 
NR  015401,  Guggenheim  Jet  Propulsion  Center,  Calif.  Inst.  Tech. 
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frequency,  l.e.  the  difference  between  the  vibrational  energies  of 
the  final  and  Initial  states,  Fy(J,K)  Is  the  rotational  energy  of 
the  Initial  state  given  by  Eq.  (8),  S(J,K)  are  the  squares  of  the 
pure  rotational  matrix  elements  given  by^ 


/  (J  +  2  ±  KHJ  -I-  1  t  K) 
2(J  +  1) 


AJ  -  +  1 


S(J,K)  = 

(  AK  =  t  1) 


/  f J  +  1  -  K)(J  1  K)(2J  +  1) 
]  2J(J  +  1) 


AJ  -  0  (15) 


( J  -  1  +  K)  ( J  4-  K) 
2J 


AJ  -  -  1 


and 


/  (J  -»■  1  +  K)(J  -<•  1  -  K) 
I  J  +  1 


S(J,K)  =  I 
(AK-0)  ^ 


(2J  +  1)K^ 
J(J  +  1) 


(J  +  K)(J  -  K) 
J 


AJ  -  +  1 

AJ  -  0  (16) 

AJ  =  -1 


and  the  amplitude  is  given  by 


I  I  '  (17) 

V  * 

where  G(v)  Is  the  vibrational  energy  of  the  Initial  state  and  R  Is 
the  purely  vibrational  matrix  element. 


O) 


N 

—  exp 
QP 


hcG(v) 


kT 
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Thus  it  follows  that  the  amplitude  Cy'  for  a  particular 
vibrational  transition  at  a  particular  temperature  is  related  to',  that 
of  any  other  vibrational  transition  at  another  temperature^  by 


C  (T)  = 

V 


V  I 


-v'* 


N(T)Q(T^) 


'*.2  (Tq)  n(T^)Q(T) 


V* 


hcG(v*) 

kT 


(18) 


It  can  be  assumed  that  N(T)  ro  (p/T)  and  the  partition  function  is 
given  by^ 


r  \ 


Q(T)  -  j  1  -  exp  i  - 


-1  /  hcB  ^  B  ^  ' 

'  ®3  I 

exp  1 - . .  !  X 


j  exp  I 


4kT 


3  ^  - 

fl(kT/hc  > 

B  B  B  ! 

®1  ®2  ®3 


,  B  ^  B  ^  ,  hcB  ^  B  '  - 
1  i  ®2  ®3  ®2  ®3 

■^12  ^  ■  ~i  kT 


exp 


.  HsGiOi 

kT 


(19) 


where  w  ^  are  the  frequencies  of  the  three  fundamental  vibrational 
modes.  This  expression  neglects  the  interaction  of  vibration  and 
rotation  and  assumes  ^Vj^^  These  assumptions  are  reasonable 

for  T  ■-  300°K. 


If  the  matrix  elements  Ry  could  be  calculated,  the  intensity 
of  any  spectral  line  could  be  deduced  from  Eq.  (14).  However,  quantum 
mechanical  expressions  for  Iff' have  never  been  derived  for  polyatomic 
molecules  such  as  H2O;  therefore  a  semi-empirical  solution  is  required. 
In  this  regard,  absolute  measurements  are  available  of  the  Integrated 
transmission  over  large  frequency  intervals  in  the  infrared.^  These 
measurements  demonstrate  that  the  vibrational  bands  of  water  vapor 


6 

J.  N.  Howard,  D.  F.  Burch  and  D.  Williams,  Sci.  Rep.  #1,  AF19(604)-516 
The  Ohio  State  University,  Dec.  1954. 
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occur  In  essentially  seven  separable  groups  In  the  frequency  range 
between  1000  and  10,000  cm"^.  Each  group  is  characterized  by  one 
or  two  relatively  Intense  vibrational  transitions  of  HHO^®  which 
dominate  all  the  Isotopic  and  higher  order  bands.  Thus  Integrated 
transmission  measurements  cannot  be  made  for  a  given  vibrational 
band  Individually,  since  such  measurements  necessarily  Include  all 
of  the  weaker  bands  In  a  given  group  vdilch  occur  In  the  same  fre¬ 
quency  Interval  with  the  most  Intense  transitions.  These  groups 
are  Identified  with  the  most  Intense  band  In  a  frequency  Interval 
and  by  the  wavelength  ^Ich  describes  the  approximate  location  of 
the  peak  absorption.  Using  the  latter  designation,  the  seven  groups 
of  the  water  vapor  spectrum  are  6.3jLi,  3.2)i,  2.7  jU.  1.87/1,  1.38/i, 
l.lju,  and  .94/i  .  The  complete  list  of  member  vibrational  transi¬ 
tions  of  these  groups  Is  given  In  Table  3. 

The  measurements  can  be  used  to  determine  one  normaliza¬ 
tion  constant  for  each  of  the  seven  groups  by  comparison  of  the 
calculated  Integrated  transmission  over  the  transitions  In  the  group 
for  a  particular  value  of  pressure  and  amount  of  water  vapor  with 
measured  values  of  this  quantity.  Details  of  this  normalization 
procedure  are  given  In  section  4.  After  this  constant  Is  determined, 
Eq.  (18)  may  be  used  to  calculate  the  Intensity  of  each  vibrational 
transition  In  the  group  provided  the  ratio  of  relative  Intensities 

(Ir^'I^/Irv'  I^)l8  known  for  the  transitions.  These  relative  In¬ 
tensities  nave  been  measured  by  several  high  resolution  studies^* 
for  many  of  the  vibrational  transitions.  For  other  bands  the  relative 
Intensities  may  be  Inferred  from  the  two  empirical  rules  discussed  by 
Herzberg:^  (1)  the  relative  Intensity  of  a  difference  band,  e.g. 

(010)  (100),  Is  approximately  equal  to  the  Boltzmann  factor  for  the 

Initial  state  times  the  relative  Intensity  of  the  corresponding 
summation  band,  e.g.  (000) —►  (110),  (2)  the  relative  Intensity  of 
bands  such  as  (010)  -*•  (110)  Is  approximately  equal  to  the  Boltzmann 
factor  for  the  Initial  state  times  the  relative  Intensity  of  corre¬ 
sponding  fundamental  bands  such  as  (000)  (100).  Furthermore  bands 

which  correspond  to  the  accidental  resonance  discussed  previously  can 
be  assumed  to  have  the  same  Intensity,  e.g.  (000)-^  (200)  and 
(000)  “♦(002).  The  relative  Intensities  of  all  the  bands  for  the 
species  HHO^^  are  given  In  Table  3  for  T  ■  300°K. 


^W.  S.  Benedict,  Bull.  Amer.  Phys.  Soc.  54  (1948). 

®W.  S.  Benedict,  Phys.  Rev.  74,  1246  (1948). 

9 

W.  S.  Benedict,  H.  H.  Claassen,  and  J.  H.  Shaw,  J.  Res.  Nat.  Bur. 
Stand.  91  (1951). 
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TABLE  3 


VIBRATIONAL  TRANSITIONS  OF  HHO^^  AT  T  -  300°K 


Group 

Transition 

Relative  Intensity 

AK 

6.3/i 

000  -*-010 

1.00 

t  1 

010  ->020 

9.6  X  10"^^^^ 

t  1 

010  100 

1.9  X  10'^^^^ 

t  1 

010  -►OOl 

-5(9) 

9.6  X  10  ^ 

0 

3.2/i 

000  -*-020 

1.00 

t  1 

010  -*-030 

4.8  X  10‘^ 

t  1 

1J\L 

000  -*-001 

1.00 

0 

000  -► 100 

.10(3) 

+  1 

010  -*011 

4.8  X  10'^ 

0 

010  -*110 

4.8  X  lO’^ 

+ 1 

1,87M 

000  -*011 

1.00 

0 

000  -* 110 

-2(7) 

2.0  X  10 

1 1 

000  -*  030 

-3(7) 

6.7  X  10 

+ 1 

010  -*021 

4.8  X  lO'^ 

0 

010  -* 120 

9.6  X  10"^ 

1 1 

010  -*040 

3.2  X  lO"^ 

1 1 

010  -*101 

4.8  X  10'^ 

0 

010  -► 200 

3.8  X  10'^ 

+ 1 

010  -*002 

5.8  X  10'3 

1 1 
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TABLE  3  (Continued) 


Group 


1.38jU 


l.lM 


Transition 

Relative  Intensity 

AK 

000-»  101 

1.00 

0 

000-^  021 

1.0  X  lO"^^^^ 

0 

000 -*•  120 

1.3  X  lO"^^^^ 

t  1 

000  -►  040 

4.3  X  10‘^ 

t  1 

000  200 

2.5  X  lO'^^^^ 

+  1 

000-*-  002 

2.5  X  10‘^ 

t  1 

010-*  111 

4.8  X  10‘^ 

0 

010 -*  031 

4.8  X  lO"^ 

0 

010-*>  130 

6.2  X  10"^ 

+ 1 

010 050 

2.1  X  10"^ 

t  1 

010-»  210 

1.2  X  10'^ 

t  1 

010-»  012 

1.2  X  10"^ 

t  1 

000  111 

1.00 

0 

000  -♦031 

-1(8) 

3.3  X  10  ^ 

0 

000  -*•  130 

-2(8) 

6.0  X  10 

1 1 

000  -♦  050 

2.0  X  10‘^ 

1 1 

000  -►210 

8.0  X  10"^^®^ 

1 1 

000  -*  012 

-1(8) 

1.2  X  10 

+ 1 

010  -►121 

4.8  X  10’^ 

0 

010  -►041 

1.6  X  lO"^ 

0 
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TABLE  3  (Continued) 


Group 


l.ljU 


.94)1 


Transition 

Relative  Intensity 

AK 

010  -♦  140 

2.9  X  10'^ 

+  1 

010  060 

9.6  X  lO'^ 

t  1 

010  -*•  220 

.5(8) 

3.8  X  10 

t  1 

010-*.  022 

.5(8) 

5.8  X  10 

t  1 

010  -*■  300 

3.2  X  10"^ 

t  1 

010 -►  003 

3.2  X  10"^ 

0 

010 -►  201 

6.4  X  10"^ 

0 

010-#-  102 

9.6  X  10"^ 

+  1 

000 121 

1.00 

0 

000  041 

.2(8) 

8.3  X  10  ^ 

0 

000 140 

1.1  X  10"^ 

+ 1 

000 -►  060 

3.7  X  lO’^ 

1 1 

000 -►  220 

2.5  X  10"^ 

1 1 

000 -♦  022 

2.5  X  10"^ 

1 1 

000 -►  300 

2.1  X  lO"^ 

1 1 

000-»  003 

2.1  X  10"^ 

0 

000  -► 102 

.1(8) 

4.2  X  10 

+ 1 

000 -*■  201 

4.2  X  10"^ 

0 

010  -►  131 

4.8  X  10"^ 

0 
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TABLE  3  (Continued) 


Group 


.94/i 


Transition  Relative  Intensity  A  K 


010  -«-051 

4.0 

X 

10"^ 

0 

010  150 

5.3 

X 

10-^ 

t  1 

010  -^070 

1.8 

X 

o 

1 

t  1 

010  -*•  230 

1.2 

X 

o 

1 

t  1 

010  -*-032 

1.2 

X 

1 

o 

t  1 

010  -*■013 

1.0 

X 

1 

o 

t  1 

010  -►  310 

1.0 

X 

1 

o 

0 

010  -►  211 

2.0 

X 

*-• 

o 

1 

0 

010  112 

2.0 

X 

io-‘ 

1 1 
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For  the  other  isotopic  species  the  relative  Intensities  are  the  values 
In  this  table  multiplied  by  the  relative  abundances. 

The  Intensity  of  any  spectral  line  can  be  obtained  from  the 
appropriate  normalization  constant,  the  data  In  Table  3  and  Eqs.  (18) 
and  (14).  However,  these  lines  correspond  to  the  equivalent  symmetric 
top  representation  of  the  water  molecule.  Thus,  It  Is  In  order  to 
determine  the  limitations  of  this  approximation  for  the  calculation 
of  the  water  vapor  transmission.  In  this  regard.  It  Is  recalled  that 
the  equivalent  S3numetrlc  top  representation  Introduces  a  double  de¬ 
generacy  Into  the  spectrum  of  rotational  energy  levels  such  that  a 
given  transition  represents  a  pair  of  transitions  In  the  actual 
asjmimetrlc  molecule  spectrum,  except  for  those  transitions  associated 
with  the  levels  K  ~  0  for  vdiich  there  Is  a  one-to-one  correspondence. 


If  the  asymmetry  Is  slight,  a  given  pair  of  transitions  will 
occur  close  together  In  frequency  and  very  near  the  frequency  of  the 
corresponding  symmetric  molecule  line.  Furthermore  the  Intensity  of 
the  symmetric  molecule  line  will  be  approximately  equal  to  the  sum  of 
the  Intensities  of  the  pair.  For  molecules  which  have  a  considerable 
as)rmmetry  such  as  H2O,  the  latter  statement  remains  valid.  This  can 
be  seen  from  the  work  of  Cross,  Halner  and  Klng.^  They  have  calculated 
and  tabulated  the  asymmetric  molecule  matrix  elements  corresponding  to 
S(J,K)  of  Eqs.  (15)  and  (16)  for  values  of  J  up  to  J  *  11.  These 
calculations  are  based  on  the  expressions  characteristic  of  an  asym¬ 
metric  molecule. 


On  the  other  hand,  the  frequencies  of  the  asymmetric  molecule 
lines  and  those  of  the  equivalent  symmetric  molecule  become  farther 
displaced  from  each  other  as  the  asymmetry  Increases.  This  Is  the  major 
difficulty  of  the  equivalent  symmetric  molecule  approximation.  Clearly 
this  displacement  of  the  spectral  lines  tends  to  define  a  lower  limit 
on  the  frequency  interval  over  which  the  approximation  may  be  used  to 
calculate  transmission.  A  determination  of  this  displacement  for  the 
lines  of  several  H2O  bands  has  been  made.  It  has  been  found  to  range 
from  one  up  to  about  100  cm“l  toward  both  lower  and  higher  frequencies. 
In  general  the  larger  values  occur  for  transitions  associated  with 
the  larger  values  of  the  quantum  number  J;  for  the  strongest  transi¬ 
tions  ( J —  4)  the  displacements  are  less  than  approximately  20  cm~^. 

Thus  It  appears  that  the  equivalent  symmetric  top  approximation  may 
certainly  be  used  to  calculate  accurate  average  transmission  values 
for  H2O  when  the  frequency  interval  of  Interest  Is  50  cm"^  or  greater. 

Fortunately  It  was  possible  in  the  present  study  to  calculate 
the  transmission  for  frequency  Intervals  much  smaller  than  50  cm"l  below 
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3400  cin~^.  Benedict,  Claassen  and  Shaw  and  Nielsen  have  measured 
and  identified  the  frequencies  and  intensities  of  the  significant 
lines  of  the  strongest  H2O  bands  corresponding  to  the  vibrational 
transitions:  (000)  —►(010)  and  (000)  —►(020)  of  the  6.3/i  and  3.2jU 
groups.  Their  calculations  were  based  on  the  asymnetrlc  top  molecule 
calculations  of  Cross,  Hainer  and  King.^  The  lines  tabulated  by 
Benedict,  Claassen  and  Shaw^  and  Nielsen^O  have  been  used  in  the 
present  study.  For  the  weaker  lines  corresponding  to  the  isotopic 
and  higher  order  bands  which  are  not  Included  in  their  tabulation, 
the  equivalent  83minetric  top  approximation  has  been  used.  Hence  in 
the  6.3/i  and  3.2/i  regions  (frequencies  below  3400  cm~^)  the  trans¬ 
missions  were  calculated  using  frequency  intervals  of  5  cm'^  (see 
Appendix  B) .  This  is  permissible  because  the  Benedict,  at  al  and 
Nielsen  tabulations  include  all  the  strong  lines,  and  the  weaker  less 
significant  lines  are  very  numerous  so  that  the  displacements  toward 
both  higher  and  lower  frequencies  tend  to  average  out. 

For  all  the  bands  occurring  beyond  3400  cm  ^  the  symnetric 
top  approximation  has  been  used.  Thus  for  frequencies  in  this  range 
the  average  transmissions  have  been  tabulated  for  20  cm~^,  50  cm~^, 
and  100  cm'l-  Intervals.  In  Appendix  A  the  coefficients  which  were 
used  in  the  quasi-random  model  have  been  tabulated  for  5  cm*^  Intervals 
for  frequencies  out  to  10,000  cm*^.  The  transitions  were  originally 
calculated  for  a  5  cm~^  interval  throughout  the  spectrum  from  this 
data.  These  results  for  5  cm~^  Intervals  were  averaged  to  obtain  the 
transmissions  over  20  cm**^,  50  crn’l-,  and  100  cm~^  given  in  appendices 
C,  D,  and  E  respectively. 

Another  difficulty  Inherent  in  the  equivalent  symnetric 
top  approximation  is  the  following:  For  the  asymnetrlc  molecule  there 
are  allowed  transitions,  but  with  low  intensity,  which  correspond  to 
the  I AK  I  —  2  transitions  for  the  symmetric  top  molecule.  These  are 
forbidden  by  the  selection  rules  for  the  symnetric  top,  Eqs.  (15) 
and  (16),  and  not  represented  vdien  the  equivalent  symmetric  top 
approximation  is  used.  However,  such  transitions  are  many  orders  of 
magnitude  less  intense  than  those  corresponding  to  |  AK|  ■  0,1.^ 

Hence,  in  general  they  contribute  insignificantly  to  the  absorption 
and  their  exclusion  is  not  expected  to  cause  serious  errors  in  the 
calculations  of  transmission. 

Finally,  it  should  be  noted  that  the  rotational  levels  of  a 
molecule  such  as  H2O  are  alternatively  symmetric  and  anti- symmetric 


^°H.  H.  Nielsen,  Phys.  Rev.  59,  565  (1941). 
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and  the  ratio  of  the  statistical  weights  of  these  two  types  of 

states  is  given  by  where  s  is  the  spin  of  the  two  identical 

8 

nuclei  in  the  molecule.  Since  the  spin  of  the  proton  is  the  anti¬ 
symmetric  levels  of  H2O  are  weighted  3  times  as  much  as  the  symmetric 
levels.  Neglecting  other  functional  dependencies,  transitions  with 
anti-s}rmmetric  initial  states  are  3  times  as  intense  as  those  with 
symmetric  initial  states.  When  the  equivalent  symmetric  top  approxi¬ 
mation  is  used,  the  intensities  given  by  Eq.  (14)  for  transitions 
with  K  or  K'  *  0  must  be  multiplied  alternatively  by  the  factors 
1  and  3,  since  there  is  a  one-to-one  correspondence  between  these  and 
the  transitions  of  the  asymmetric  molecule.  When  K  and  K'  are  greater 
than  zero,  a  given  equivalent  symmetric  molecule  transition  corresponds 
to  a  pair  of  transitions,  one  symmetric  and  the  other  anti-symmetric, 
of  the  asymmetric  molecule  and  thus  is  weighted  by  the  factor  4. 

For  the  HDO  molecule,  this  intensity  alternation  does  not  occur  and 
the  weighting  factors  are  taken  to  be  2,  2,  and  4  rather  than  1,  3, 
and  4. 


I.  Schiff,  Quantum  Mechanics,  McGraw-Hill,  New  York  (1949). 
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SECTION  3 


THE  QUASI- RANDOM  MODEL  OF  BAND  ABSORPTION 


3.1  THE  QUASI- RANDOM  MODEL 

Several  different  spectral  models  have  been  derived  which 
represent  the  average  transmission  of  radiation  absorbed  by  a  group 
of  spectral  lines  which  occur  within  a  finite  frequency  interval. 
However,  the  previous  models  often  do  not  adequately  represent  the 
actual  transmission  over  a  wide  range  Of  pressures  and  path  lengths. 
The  first  difficulty  is  that  the  assumed  arrangement  of  the  spectral 
lines  within  the  frequency  interval  does  not  correspond  adequately 
with  the  actual  spectrum.  Secondly,  the  contributions  from  lines 
lying  outside  the  interval  are  incorrectly  taken  into  account. 
Furthermore  the  variation  in  Intensity  of  the  lines  is  usually 
represented  by  Inaccurate  approximations.  Other  approximations 
which  are  only  valid  in  certain  limited  cases  are  sometimes  intro¬ 
duced  in  order  to  simplify  the  calculations. 

12  3 

For  example,  the  statistical  model  ’’  is  based  on  the 
assumption  of  a  random  arrangement  of  the  spectral  lines  fdiile  the 
Elsasser  model  assumes  a  constant  spacing  between  the  lines.  For 
the  statistical  model  an  expression  for  the  transmission  is  readily 
derived  from  the  assumption  that  each  line  center  has  an  equal 
probability  of  falling  at  any  frequency  within  the  interval.  Itow- 
ever,  in  the  spectrum  of  an  actual  molecule  a  given  frequency 


^H.  Mayer,  Methods  of  Opacity  Calculations,  Los  Alamos,  LA-647, 
Oct.  31,  1947. 

^R.  M.  Goody,  Quart.  J.  Roy.  Met.  Soc.  T^,  165  (1952). 

^G.  N.  Plass,  J.  Opt.  Soc  Amer.  690  (1958);  W,  868  (1960). 
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Interval  is  generally  characterized  by  the  occurrence  of  a  relatively 
small  number  of  stronger  lines.  These  strong  lines  occur  neither  at 
random  nor  are  they  regularly  arranged.  A  model  Is  required  which 
localizes  the  stronger  spectral  lines  so  that  their  positions  rela¬ 
tive  to  each  other  and  the  Interval  are  taken  Into  account.  The 
quasi-random  model  localizes  the  lines  to  within  a  frequency  Interval 
\^lch  Is  small  compared  to  the  final  Interval  of  Interest. 

The  quasi-random  model  Is  characterized  by  the  following 
six  features  which  are  discussed  In  more  detail  later  In  the  report. 

1.  The  stronger  spectral  lines  are  well  localized  by  the 
use  of  small  frequency  Intervals  of  width  D.  Within  this  small 
Interval  D  the  lines  are  assumed  to  be  arranged  at  random.  The 
transmission  over  larger  frequency  Intervals  A  v>  >D  Is  calculated 
by  averaging  the  results  from  the  smaller  Intervals. 

2.  The  removal  of  artificial  effects  due  to  the  arbitrary 
division  of  the  spectrum  Into  frequency  Intervals.  This  Is  done  by 
the  displacement  In  frequency  of  the  mesh  which  divides  the  spectrum. 
The  transmission  Is  first  calculated  for  one  mesh,  which  defines  a 
certain  set  of  frequency  Intervals.  The  calculation  Is  repeated  for 
another  mesh  vdilch  Is  shifted  In  frequency  by  %D  relative  to  the  first 
set.  In  principle,  the  mesh  can  be  shifted  n  times  by  an  amount  D/n. 
The  final  calculated  transmission  Is  the  average  of  the  results  for 
each  of  these  meshes. 

3.  The  transmission  Is  calculated  from  expressions  which 
are  valid  for  a  finite  and  possibly  small  number  of  spectral  lines 
In  the  Interval.  The  expression  is  not  used  which  Is  valid  In  the 
limit  as  the  number  of  lines  becomes  very  large. 

4.  The  spectral  lines  In  each  frequency  Interval  are  divided 
Into  Intensity  subgroups  vdilch  are  fine  enough  to  simulate  the  actual 
Intensity  distribution.  In  the  calculation  the  average  Intensity  of 
all  the  lines  In  each  subgroup  Is  used  together  with  the  actual  number 
of  lines. 


4 

5.  The  Landenberg-Relche  approximation  assumes  that  the 
frequency  Interval  Is  infinitely  wide.  This  expression  cannot  be  used 


Landenberg  and  F. 


Reiche,  Ann.  Physlk  181  (1913). 
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In  this  calculation  becauae  of  the  small  frequency  Intervals  into 
which  the  spectrum  Is  divided.  Instead  more  general  expressions 
have  been  derived  which  are  valid  for  a  finite  frequency  Interval, 
and  these  are  used  with  the  quasi- random  model. 

6.  The  contribution  from  the  wings  of  the  spectral  lines 
whose  centers  are  outside  a  given  Interval  Is  Included.  The  wing 
effect  Is  treated  In  the  same  detailed  manner  as  the  contribution 
from  the  lines  within  the  Interval. 

The  grouping  of  spectral  lines  within  small  frequency 
Intervals  Is  basic  to  the  quasi-random  model.  If  these  Intervals 
are  chosen  small  enough,  a  transmission  characteristic  of  each 
Interval  can  be  obtained  by  assuming  a  random  position  for  the  lines 
within  the  small  Interval.  An  Integrated  average  transmission  over 
much  larger  frequency  Intervals  can  be  obtained  by  averaging  the 
results  over  the  small  Intervals. 

As  the  size  of  a  frequency  Interval  Is  decreased,  the 
relative  position  of  the  lines  with  respect  to  one  another  within 
the  Interval  becomes  less  and  less  Important.  Each  line  may  be 
said  to  be  localized  within  an  error  defined  by  the  Interval  size,  D. 
The  lines  within  the  small  Interval  D  are  considered  to  occupyrandom 
positions  within  the  Interval.  This  assumption  cannot  be  made  for 
larger  frequency  Intervals  Av  idiere  the  random  hypothesis  can 
Introduce  errors  for  several  reasons.  First,  the  spacing  between 
the  lines  may  be  fairly  uniform.  In  which  case  the  random  hypothesis 
overestimates  the  overlapping  of  the  lines.  Second,  the  strong  lines 
may  all  occur  near  one  end  of  the  larger  frequency  Interval.  In  this 
case  for  large  path  lengths  the  random  hypothesis  might  predict 
almost  complete  absorption,  whereas  In  actuality  only  a  small  part 
of  the  Interval  In  the  region  of  the  strong  lines  Is  strongly  ab¬ 
sorbing  and  the  rest  of  the  Interval  Is  nearly  transparent. 

In  both  of  these  cases.  If  the  original  Interval  Is  sub¬ 
divided  Into  sufficiently  small  Intervals,  each  of  which  may  contain 
only  one  or  at  most  a  few  strong  lines,  the  transmission  can  then 
be  calculated  accurately  by  assuming  the  random  placement  of  the 
spectral  lines  within  the  small  Interval.  Of  course.  If  the  lines 
are  arranged  nearly  at  random  In  the  original  large  Interval,  then 
nothing  Is  lost  by  a  subdivision  Into  finer  Intervals.  Since  the 
model  localizes  the  lines  within  the  smaller  frequency  Interval  D, 
while  calculating  the  absorption  for  a  larger  Interval  Av  ,  It  Is 
called  quasi-random. 
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If  the  frequency  spectrum  Is  divided  Into  small  Intervals, 
D,  the  transmission  T  in  each  of  these  small  Intervals  is  calculated 
from  the  equations  of  the  statistical  model. If  there  are  nj^ 
lines  in  the  Interval  D|(  with  their  line  centers  at  the  frequencies 

(i  «  1,2 . n)() ,  then  the  transmission  at  the  frequency  v  is 

affected  by  these  nj^  lines  so  that 


(V) 


“k 


-S^ub(  V, 
e 


^i> 


d  V 


(20) 


where  is  the  intensity  of  the  1^  line,  u  is  the  amount  of  absorb¬ 
ing  gas  per  unit  area,  and  b  is  the  line  shape  factor.  The  trans¬ 
mission  is  calculated  at  some  frequency  v  vdilch  for  convenience  is 
usually  taken  at  the  center  of  the  interval  Since  this  trans¬ 

mission  corresponds  to  the  average  of  all  permutations  of  the  posi¬ 
tions  of  the  spectral  lines  within  the  interval,  it  is  also  assumed 
to  represent  the  most  probable  transmission  for  the  entire  interval. 

In  principle  the  transmission  of  each  spectral  line  should 
be  calculated  from  its  intensity  and  position  and  then  substituted 
into  Eq.  (20).  Because  of  the  very  large  number  of  spectral  lines 
which  were  considered  in  the  present  calculation,  it  was  not  practical 
to  do  this  evjen  with  a  presently  available  high  speed  digital  conqputer. 
Thus  some  method  must  be  found  to  simulate  the  actual  intensity  dls-  ’ 
tribution.  A  coiiq>arison  of  calculated  intensity  distributions  with 
those  >dilch  have  been  suggested  for  use  with  the  statistical  model^>^ 
shows  that  there  is  usually  poor  agreement  between  the  two. 

The  procedure  for  intensity  simulation  adopted  here  is  to 
divide  the  lines  in  each  frequency  interval  into  subgroups  by  intensity 
decades.  Vftien  the  line  intensities  are  calculated,  the  results  are 
grouped  according  to  these  intensity  decades.  The  average  intensity 
of  the  lines  in  each  decade  is  then  calculated. 

Test  calculations  have  shown  that  only  the  first  five 
intensity  decades  for  each  frequency  interval  influence  the  value  of 
the  transmission.  Thus  only  the  data  for  the  five  strongest  intensity 
decades  need  be  retained  for  each  frequency  interval.  It  should  be 
emphasized  that  the  numerical  range  of  the  intensity  decades  retained 
are  quite  different  for  a  frequency  interval  which  contains  strong 
lines  than  for  one  which  has  weak  lines.  The  criterion  is  applied 
separately  to  each  frequency  interval. 
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Thus  the  transmission  can  be  calculated  from  the  average 
value  of  the  intensity  together  with  the  number  of  lines  n^  in 
each  intensity  decade  from  the  equation 


T  (V)  -  TT 


1 

D  a 


pfc  If 


(21) 


The  number  of  lines  n^^  in  the  frequency  interval  is  given  by 


5 

n.  =  E  n.  (22) 

k  i 


Of  course,  the  lines  could  be  grouped  by  Intensity  using  a  smaller 
or  larger  factor  than  the  decade  of  intensity  adopted  here.  However, 
this  choice  seemed  reasonable  since  it  yields  intensity  subgroups 
idilch  are  sufficiently  small  Intervals  so  that  the  average  value  is 
never  much  different  from  the  actual  intensity  of  the  lines  in  the 
decade,  while  on  the  other  hand  the  interval  is  sufficiently  large 
so  that  the  number  of  different  variables  used  in  the  calculation 
are  still  within  the  memory  limitations  of  a  modern  digital  computer. 
All  of  the  lines  v^lch  occur  in  a  particular  intensity  group  are 
considered  to  occur  at  random  in  the  frequency  interval  Dj^,  as  is 
indicated  by  Eq.  (21).  Furthermore  each  intensity  group  is  treated 
Independently  of  any  other. 

The  integral  in  Eq.  (21)  represents  the  transmission  over 
the  finite  interval  D)^  of  a  single  spectral  line  representative  of  the 
average  intensity  of  the  decade.  For  the  small  intervals  used  in  the 
present  calculation,  it  is  important  to  separate  the  absorption 
caused  by  a  given  spectral  line  into  two  parts:  the  absorption  vdilch 
it  produces  within  the  interval  which  contains  its  line  center  and  the 
absorption  produced  by  it  in  all  other  frequency  intervals  arising 
from  its  wings.  For  the  calculation  of  the  absorption  within  the 
frequency  interval  tdiich  contains  the  line  center  the  Landenberg- 
Reiche^  equation  cannot  be  used.  It  assumes  an  infinite  frequency 
interval,  whereas  the  intervals  in  the  present  calculation  are  not 
necessarily  large  compared  to  the  effective  width  of  the  line. 

Although  equations  for  the  absorption  of  lines  over  finite  Intervals 
have  been  obtained  previously  in  certain  limits,  a  general  expression 
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has  not  been  available.  !nils  result  Is  given  In  section  3.2,  Eqs. 
(27),  (30),  and  (32).  Because  of  the  small  number  of  strong  lines 
in  the  smaller  frequency  Intervals,  It  Is  Important  to  calculate 
the  transmission  of  each  line  which  represents  an  Intensity  decade 
and  then  to  multiply  these  transmissions  as  Indicated  In  Eq.  (21). 
Accordingly,  the  exponential  expression  for  the  transmission  de¬ 
rived  from  the  statistical  model  tdien  the  number  of  lines  approaches 
Infinity  may  not  be  used  since  the  number  of  lines  In  an  Intensity 
decade  In  these  smaller  frequency  Intervals  may  be  small. 

The  transmission  at  a  frequency  v  Is  also  Influenced  by 
the  wings  of  spectral  lines  whose  centers  are  outside  of  the  Interval 
vdilch  contains  the  frequency  v  .  The  wing  transmission  Tni  from 
each  of  these  Intervals  Di  Is  calculated  by  assuming  a  random  dis¬ 
tribution  of  the  lines  within  the  Intervals  D^.  Thus  by  the  random 
hypothesis  the  transmission  at  the  frequency  v  Is  the  product  of 
all  transmissions  calculated  for  the  lines  with  centers  In  other 
Intervals  taken  with  the  transmission  for  the  given  Interval  Itself, 
so  that 


T(  V) 


00 

Tt  T  (v) 
n. 


(23) 


where  the  subscript  n^  Indicates  that  there  are  nj^  lines  In  the 
frequency  Interval  D^.  The  transmission  for  the  lines  whose  centers 
are  In  the  Interval  are  given  by  Eq.  (21).  The  lines  In  suffi¬ 
ciently  distant  Intervals  from  the  particular  frequency  v  have  a 
negligible  effect  on  the  transmission,  so  that  In  practice  the 
product  (Eq.  23)  usually  needs  to  be  evaluated  only  over  a  few  terms. 
Approximate  expressions  for  the  wing  transmission  Tni  are  given  In 
section  3.2,  Eqs.  (27)  and  (29)  for  the  Lorentz  line  shape  and  In 
section  3.3,  Eq.  (42),  for  the  Benedict  modification  of  the  Lorentz 
line  shape. 

Finally,  It  Is  necessary  to  remove  the  artificial  effects 
arising  from  the  arbitrary  division  of  the  spectrum  Into  frequency 
Intervals.  In  order  to  Illustrate  this  effect  consider  a  set  of 
six  equally  Intense  lines  spanning  a  frequency  range  of  10  cm~^. 

If  the  lines,  for  simplicity,  be  nearly  evenly  spaced,  the  trans¬ 
mission  calculated  at  the  center  of  the  two  Intervals  containing 
them  might  be  expected  to  be  equal.  This  would  be  the  case  If  three 
lines  fell  Into  each  5  cm'^  interval,  for  the  model  subsequently 
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attributes  a  random  distribution  to  each  set  of  lines  falling  within 
the  chosen  intervals.  But  since  the  uniform  mesh  is  started  at  some 
arbitrary  position,  it  is  conceivable  that  two  lines  might  fall  into 
one  interval  and  four  into  the  other.  This  situation  is  illustrated 
in  Fig.  3.1  where  the  mesh  labeled  "A"  represents  the  former  case 
while  "B"  Illustrates  the  latter.  The  six  bars  schematically  repre¬ 
sent  the  lines  and  the  arrows  the  frequency  divisions. 

The  transmissions  calculated  at  the  centers  of  the  Intervals 
in  case  "B"  differ  appreciably  from  those  in  case  "A",  whereas 
physically  one  would  expect  them  to  be  nearly  equal  for  lines  with 
such  uniform  spacing.  The  disparity,  of  course,  may  be  blamed  firmly 
upon  the  positioning  of  line  Q)  .  As  far  as  the  model  is  concerned, 
the  fact  that  line  Q)  lies  very  close  to  the  artificial  mesh  boundary 
is  of  no  consequence;  in  case  "A"  it  has  a  random  probability  of  being 
anywhere  within  interval  I  whereas  in  case  "B"  it  may  lie  anywhere 
within  II'.  Thus  the  model  does  not  take  into  account  such  borderline 
situations.  The  Introduction  of  such  an  error  is  of  no  practical 
importance  if  the  results  of  a  number  of  frequency  intervals  are  to 
be  averaged  to  obtain  the  integrated  average  transmission  over  a  larger 
Interval.  However,  if  the  detailed  average  effective  transmission  is 
required  over  smaller  frequency  intervals,  then  the  model  must  be 
suitably  modified. 

The  modification  of  the  transmission  calculations  necessary 
to  remedy  these  artificial  errors  is  extremely  sinqple.  The  transmission 
is  calculated  twice  using  tw>  meshes  one  of  %^lch  has  been  displaced 
D/2  cm~^  with  respect  to  the  other.  Thus  the  first  calculation  yields 
values  at  v^,  +  D,  +  2D,  ...  while  the  second  produces  trans¬ 

mission  values  at  Vq  +  D/2,  Vq  +  3D/2,  Vq  +  5D/2,  ...,  where  Vq 
corresponds  to  the  lowest  frequency  at  vdiich  the  transmissions  are 
calculated.  Thus  transmission  values  are  calculated  at  the  following 
frequencies: 


V 


1 


V 


2 


V  .  V  _  +  D 

J  o 

v^  +  3D/2 


Vj  V^  +  2D 
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The  correct  average  transmissions  at  a  particular  frequency  sub¬ 
sequently  are  calculated  from  these  values  by  the  following  simple 
formula: 


T( 


T(v^-l)  +  T(v^)  + 
3 


(24) 


where  the  T's  refer  to  the  effective  transmissions  calculated  by  the 
dual  mesh  scheme  and  labeled  as  shown  In  the  table  above.  Such  a 
procedure  automatically  compensates  for  the  artificial  errors  Intro¬ 
duced  by  the  arbitrary  mesh  positioning  and  provides  an  excellent 
representation  of  the  average  transmission  over  small  frequency  Inter¬ 
vals.  All  of  the  average  effective  transmission  values  given  In  this 
report  are  averaged  In  this  manner  using  two  meshes  shifted  relative 
to  each  other  by  ^D. 


3.2  THE  TRANSMISSION  FOR  A  FINITE  FREQUENCY  INTERVAL 

In  this  section  an  expression  Is  derived  for  the  transmission 
of  a  finite  frequency  Interval  for  spectral  lines  with  the  Lorentz  line 
shape.  A  finite  number  of  spectral  lines  are  assumed  to  be  In  the 
Interval.  Expressions  have  been  derived  previously  for  the  transmission 
according  to  the  statistical  model  for  a  finite  number  of  spectral 
llnes^*2,3  and  In  certain  limits  for  a  finite  frequency  Interval.^  How¬ 
ever,  In  practical  calculations  which  use  the  statistical  model  tt^ 
assumption  seems  to  have  always  been  made  that  the  equations  for  an 
infinite  number  of  spectral  lines  and  an  Infinite  frequency  Interval 
are  sufficiently  accurate.  Such  expressions  cannot  be  used  for  the 
small  frequency  Intervals  which  are  used  In  the  present  calculation. 

Therefore  a  general  expression  for  the  transmission  for 
finite  frequency  Intervals  and  a  finite  number  of  spectral  lines  Is 
derived.  In  this  section  It  Is  assumed  that  the  Lorentz  line  shape 


b(  V  V.) 
»  * 


a/7T 

(V  -  + 


(25) 


describes  the  actual  line  contours,  where  a  Is  the  pressure  broadened 
half-width  of  the  line  whose  center  lies  atv^.  VRien  pressure  broaden¬ 
ing  Is  the  dominant  effect,  the  lines  of  H2O  are  well  described  by 
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the  Lorentz  line  shape.  For  CO2  the  wings  of  the  lines  seem  to  be 
described  better  by  the  Benedict  modification  of  the  Lorentz  line 
shape. 5  This  is  described  further  in  section  3.3.  Consider  a  fre¬ 
quency  interval  D  which  lies  between  the  frequencies  Vg  and  Vg  +  D. 
The  following  symbols  are  used  in  the  derivation  of  the  transmission 
over  a  finite  frequency  interval: 


p  ="  2a  /D 


y  -  -  D/2 

z  =  V  -  V  -  D/2 

o 

T)  -=  2y/D 

e  *  2z/D 


(26) 


Thus  for  the  Lorentz  line  shape  the  transmission,  T  ,  is  given  by 
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(27) 


Define  by  the  expression 


"  ^1,  “/P 


S,  p  u 
k'^o 

TTtt  p 
O*^ 


(28) 


where  a^  is  the  half-width  at  the  pressure  p^^. 


At  a  particular 


5 


W.  S.  Benedict,  private  communication. 
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temperature,  the  coefficients  ^  need  be  calculated  only  once  for 
each  frequency  Interval  and  Intensity  group.  The  coefficients 
^  together  with  the  number  of  spectral  lines  for  each  of  the  top 
five  intensity  groups  represent  sufficient  information  about  the 
spectrum  for  the  calculation  of  the  absorption  in  a  frequency  interval. 
When  ^  is  known,  then  the  coefficients  for  any  pressure  and  path 
length  combination  may  be  readily  calculated  from  Eq.  (28).  Indeed, 
this  is  one  of  the  most  important  simplifications  of  the  present  model. 

Because  p  =  2a/D  is  usually  very  small,  Eq.  (23)  can  be 
divided  easily  into  two  distinct  parts.  There  is  the  direct  contri¬ 
bution  from  the  lines  whose  centers  are  within  the  given  interval  in 
addition  to  the  wing  contribution  from  the  lines  in  all  other  intervals. 


So  long  as  V  is  not  chosen  too  close  to  the  edge  of  an 
interval,  which  is  generally  taken  to  be  large  compared  to  the  half- 
width  of  any  line,  then  (£-'n)^»P^  for  all  the  factors  in  Eq.  (23) 
which  correspond  to  lines  whose  centers  are  in  other  frequency 
intervals.  Hence  for  the  wing  contributions  p^  may  be  neglected 
compared  to  (e-r))2.  Thus  the  integral  in  Eq.  (27)  can  be  written  as 
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.  V'  A  ,  '/a 

erf  —  -  erf  ^ 


where  A  =  p  and  €  >1.  Thus  Eq.  (29)  represents  a  readily  evaluated 
quantity. 


The  direct  contribution,  on  the  other  hand,  is  more  difficult 
to  evaluate.  Let 
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vdiere  6  =  p/a.  Therefore  returning  to  the  direct  contribution  one 
finds 
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In  section  3.4  it  is  shown  that  for  the  general  case 


Q(  ?,6)  -  ^  e'  [lo(?/2)  +  l^(!i/2) 


iTT  -f) 


t  00 

+  26e'^^^  ^  I^(^/2)8lnn^ 

n*l 


(32) 


+  6?  ^ 

n«l 


sin 

n 


where  f  ^  2  tan  (>  ,  5  »  p/a  and  the  I  's  are  the  Bessel  functions 
of  imaginary  argument.  Although  the  sums  in  Eq.  (32)  are  rapidly 
convergent,  in  the  two  cases  when  ^  »  1  and  ^  <  1  useful  alter¬ 

native  forms  may  be  found,  whose  derivation  is  also  found  in  section  3.4. 
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Specifically  for  |  <  1  one  has 


0(^.6  ) 


1 


00 


(-1) 


1 


n+1 


(n+l)! 


J 


where 


6 

2n(l+6^)” 


+ 


2n-l 

2n 


-  tan‘\l/6  ) 


(33) 


(34) 


The  series  In  Eq.  (34)  Is  absolutely  convergent  for  all  4  but  Is  more 
useful  when  4  Is  small. 

When  ^  Is  large,  on  the  other  hand,  an  asymptotic  repre¬ 
sentation  Is 


P(^,  6)  -  e*^^  [l^,(l/2)  +  I^(^/2) 


+  2 


m  m 


m“0 


m. 


+  2  6  e  R 


(35) 


where 


'=^m  ^  2|  ^m-1 


'4)  “  ^ 


(36) 


e  -  6  /V/l+6 
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with  the  remainder  term,  Rp,  given  by 


Jy2Pe-^y\y  r 

-U  Jo  (l-y^)P“’ 


(37) 


and  the  Barnes'  coefficient  (  Ol)^  defined  by 


(a)jj  ■  a(  a+1)*  •  •  (  a+n-1),  (  a)jj  ■  1 


(38) 


On  the  basis  of  Eqs.  (29)  and  (31)  the  effective  trans¬ 
mission  at  a  particular  frequency  may  now  be  written  as 


T(v)  -  3:  TT  \ 

^  ll‘k  ^ 

where  the  subscript  k  refers  to  the  particular  Interval  In  which 
lies  and  Ti((  V)  Is  obtained  from  Eqs.  (21),  (27),  (31),  and  (32) 
while  Ti(v),  for  1  k.  Is  readily  calculated  from  Eqs.  (21),  (27) 
and  (29). 

The  water  vapor  and  carbon  dioxide  (SSD-TDR-62-127  *■  Volume  III) 
transmissions  were  all  calculated  on  the  basis  of  this  formalism.  The 
absorption  from  a  single  spectral  line  which  represents  an  Intensity 
group  In  a  finite  frequency  Interval  was  calculated  from  these  equations. 
The  other  lines  In  a  small  frequency  Interval  are  assumed  to  be  dis¬ 
tributed  at  random  so  that  the  transmission  of  the  various  lines  can 
be  BHiltlplled  together  according  to  Eq.  (21).  The  transmissions  cal¬ 
culated  from  the  wing  absorption  of  lines  outside  the  Interval  are 
multiplied  according  to  Eq.  (39).  Finally  the  results  are  calculated 
for  two  different  frequency  meshes  shifted  an  amount  with  respect 
to  each  other  and  these  results  are  averaged  to  obtain  the  average 
effective  transmission  for  the  Interval. 

3.3  BENEDICT  MODIFICATION  OT  THE  LORENTZ  LINE  SHAPE 

As  mentioned  In  section  3.2,  Benedict^  has  proposed  a  modifi¬ 
cation  of  the  Lorentz  line  shape  for  OO2  on  the  basis  of  his  experimental 
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examination  of  the  absorption  from  some  particular  lines.  He  has 
found  that  the  line  shape  in  the  wings  drops  off  much  faster  than 

1^  '  ^ol~^  required  by  the  Lorentz  line  shape.  On  the  basis  of 
these  findings  he  has  suggested  the  following  seml-emplrlcal  line 
shape 


b(  V  ,  vp 


g  /  7T 


2  2 
( V- v^y+ a 


,  I V-  V  I  <  d 


(40) 


A  a/7T  exp-a|  v  - 

2 

(.V-  V  +  a 


,1  V  -  V  ^1  >  d 


where  a  ■  0.0675,  b  *  0.7,  d  »  5.0  cm 
makes  the  two  forms  continuous  at  |  V  -  V 


and  A  Is  a  constant  which 
^  I  *  d  “  5 . 0  cm"  1 . 


The  effective  transmission  is  now  calculated  In  a  manner 
similar  to  that  described  In  section  3.2  with  one  Important  change. 
Since  the  frequency  Interval  chosen  for  these  calculations  Is  5  cm'^, 
the  Lorentz  line  shape  Is  still  valid  for  lines  whose  centers  are 
within  the  Interval.  This  contribution  can  still  be  calculated  as 
before  from  Eqs.  (21),  (27),  (31),  and  (32).  On  the  other  hand, 
the  wing  contributions  must  be  treated  differently.  Introducing  the 
symbols  defined  In  Eq.  (26),  the  contributions  from  the  wings  will 
involve  terms 


'w  <■')  --B  J 


exp 


ASj^ug  exp-a  I V  - 

7rf(v  -  V  )^+  g^J 


lb  1 


dv 


2  4 


exp  - 


A^  p^exp-  [a(D/2)'’le  -  T)  | j 


(41) 


(e  -  7i)' 


dT) 


2  2 
where  p  has  again  been  neglected  compared  with  ( e -  7) )  (see 

section  3.2).  Equation  (41)  cannot  be  reduced  to  an  analytic  form 

as  was  done  for  the  transmission  from  the  wings  of  a  line  with  the 
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Lorentz  line  shape,  Eq.  (29).  Since  the  absorption  drops  off  even 
more  rapidly  with  frequency  in  the  Benedict  modification  than  with 
the  Lorentz  line  shape,  the  wing  contributions  may  be  calculated 
to  a  high  degree  of  accuracy  by  assuming  that  the  lines  are  all 
located  at  an  average  position  rather  than  a  random  one.  The  average 
position  chosen  was  the  center  of  the  interval,  i.e.  7)  -  0.  In 
this  event  Eq.  (41)  becomes 


T  (V) 


exp 


.} 


A^jjP^exp-  [a(D/2)'*  I  e|  **] 


(42) 


The  effective  transmission  for  the  Benedict  modification  of  the 
Lorentz  line  shape  at  a  particular  frequency  V  is  now  written  as 


T(V  )  -  T.(v)  7[  T.(v) 

iA  ^ 


(43) 


where  the  subscript  k  refers  to  the  interval  in  which  v  lies, 

T|c(v)  is  calculated  from  Eq.  (31)  and  the  Ti(v)'s  are  calculated 
on  the  basis  of  Eq.  (42). 

The  CO2  transmissions  were  all  calculated  on  this  basis. 
Nevertheless,  for  comparison  purposes,  they  were  recalculated  using 
the  Lorentz  line  shape  at  300°K  and  one  atmosphere  pressure.  Appendix 
D  of  SSD-TDR-62-127  -  Volume  III  gives  the  results  of  this  calculation. 
It  is  interesting  to  note  that  there  is  a  very  small  difference 
between  the  two  calculations  even  for  large  path  lengths. 


3.4  DERIVATION  OF  THE  TRANSMISSION  EQUATIONS  FOR  THE  DIRECT 
CONTRIBUTIONS 


Consider  now  the  evaluation  of  the  integral  occurring  in 
Eq.  (30),  viz. 


.6) 


dy 


(44) 
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1 


Letting  y 


6^/T-^/z,  Eq. 


(44)  becomes 


Q(|.  6) 


(45) 


Integrating  by  parts  once 

Q  (  ^  .  5  )  -  e'^  ^  ^-26^  /Vl-z2  e' ^  ^  dz  (46) 

6  /  Vl+  62 


I 


Consider  now 


(47) 


G 


3 


Letting  z  >  sin  0,  the  first  Integral  becomes 

G. 


I 


cos  9  e  d9 


■i  n 


I 


TT 


cos  ^ 

(1+cos  p)  e  dp 


TTe 


-^/2 


O)  -  IJj(^) 


(48) 
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where  the  substitutions  ^  ■  29  and  j8  “  have  been  made  to  obtain 
the  final  result  involving  the  usual  Bessel  functions. 

Re-expressing  the  result  in  terms  of  Bessel  functions  of 
imaginary  argument 


r  - 

■  4 


■in 


I^(?/2)  +  I^(|/2) 


(49) 


Letting  e  =  sin  ''6  /Vi+62 


-1 


tan’^  6,  the  second  integral 


in  Eq.  (47)  becomes 


■r 


2  _  -  I*  sin  0 

cos  0  e  ’  d0 


(50) 


where  again  the  substitution  z  -  sin  0  has  been  made.  Once  more  let 
d  =  2  0  to  obtain 


G  =  ^  e- 
4  ® 


L 


2e 


,  , . , ^  £cos  ^  , 

(cos  d+1)  e  ^ -  d^ 


Letting  ^  J 


2e 


i  cos  i 


d^ 


one  has 


G-i  =  T  e'  ■  G,  + 

3  4 


dG, 


I., 


(51) 


(52) 


and  hence  it  remains  only  to  evaluate  G^.  Substituting  z  >  i^/2, 
one  has 


L 

i 


-iz  cos 


2e 

cos(z  cos 


^)di  - 


i 


2e 

sin(z  cos  ^)d^ 


(53) 
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However  co8(z  cos  and  sln^z  cos  i)  may  be  represented  by  their 
absolutely  convergent  Jacobi”  expansions,  viz. 

« 

cos(z  cos  ^)  ■  J  (z)  +  2  J]  (-1)”  J«  (z)  cos  2n  6 
o  ",  zn 

n-1 

00  (54) 

sin(z  cos  (i)  ■  2  ^  (-I)*'  ‘^2n+l^*^  cos(2n+l)  ^ 

n«o 


Hence 


J^(z)  i/+  2 


I 


(-1)“  J2n(z) 


sin  2n  il/ 
2n 


(.1)^  J  sln(2nf^V^ 

^  ^  '^2n+r  ^  2n+l 


(55) 


where  if/ 


2e 


2  tan  ^  6 


Re-expresslng  In  terms  of  ^  one  has 


Gj  -  ye  /2)  *  +  2  f  /2)  liEJai. 

n*l 


w 

+  2  ^  I-  ^  ( 


8ln(2n-H)  f 


n»o 


2n+l'^^^^  2n+l 


(56) 


OO 

1(^/2)  +  2  Z  1  (  i  12)  S.yFi^J.. 


n«l 


^C.G.J.  Jacobi,  Journal  fur  Math.  15,  12  (1836). 
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Therefore 


( 


j  /2)  ^  +  Ij(  e/2)  Tjj 

[ 

+  -|-  ^  I^(e/2)  sin 
n*l 
00 

+  2^  +  l^^^(e/2)] 


1  .-4/2 
'3  4 


00 

4  V" 


(57) 


If  Eqs.  (49)  and  (57)  are 
previously  quoted  result. 


substituted  into  Eq.  (46)  and  (47),  the 
Eq.  (32),  is  obtained. 


Although  Eq.  (32)  is  absolutely  convergent,  vrtien  ^  <  1 
and  4  >>  1  two  simplified  forms  may  be  obtained  (Eqs.  33  and  35) 

!f%^nnr  Eq-  W  directly  by  expanding  the  kernel 

integrating  term  by  term.  This  procedure 
immediately  produces  Eq.  (33)  which  is  itself  absolutely  convergent. 


When  ^  »  1,  an  asymptotic  form 
from  Eqs.  (46)  and  (47)  and  expanding  Vl  - 
via  the  binomial  expansion  theorem.  Thus 


nifiy  be  obtained  by  starting 
22  under  the  integral  sign 


Q(4,6)  -  e  )  .  ^  e"  [i^(^/2)  +  l^(^/2)] 

6  /  V 1+  6^ 


+  2  6^ 


•a2  e'^"  dz 


(58) 


But 


P-1 
dz  -  1C 


m*o 


_  m  m 


+  R 
P 


(59) 


3-20 


where<<J  is  given  by  Eq.  (36)  and  is  the  remainder  term  given  by 
Eq.  (37).  The  summation  represe/its'^the  first  p  integrals  obtained  by 

2p-2 


expanding  Yl-z  to  order  z  The  remainder  term  was  obtained 

on  the  basis  of  a  formalism  apparently  due  to  Graf  and  Gubler.^ 


The  second  term  in  Eq.  (58)  is  the  familiar  Landenberg 
and  Reiche  result^  which  is  valid  idien  the  exponential  and  Integral 
term  in  this  equation  can  be  neglected.  The  expression  for  the  absorp¬ 
tion  over  a  finite  interval  \dien  the  strong  line  approximation  is 
valid^  is  obtained  from  Eq.  (58)  by  inserting  an  asymptotic  expression 
for  the  Bessel  functions  and  evaluating  the  leading  term  of  the  inte¬ 
gral  from  Eqs.  (59).  and  (36).  The  absorption  over  a  finite  Interval 
vdien  the  weak  line  approximation  is  valid^  is  obtained  from  Eqs.  (33) 
and  (34) . 


^J.  H.  Graf  and  E.  Gubler,  Einleitung  in  die  Theorie  der  Besselshen 
Funktlonen,  86  (Bern,  1898). 
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SECTION  4 


SPECIFIC  DETAILS  OF  THE  CALCULATION 


Prior  to  calculating  the  transmissions  based  on  the  pre¬ 
viously  described  model,  a  coefficient  table  was  generated  for  the 
two  frequency  meshes  at  each  of  the  temperatures  considered.  Each 
of  these  coefficient  tables  ia^composed  of  ^  ^  function 

of  frequency  Interval,  where  nj^  is  the  number  of  lines  in  the 
intensity  decade  and  is  defined  by  Eq.  (28).  At  each  frequency 

Interval  the  coefficients  were  calculated  for  up  to  five  intensity 
decade  groups  in  the  following  manner: 

(1)  The  vibrational  intensity  factors  C^'of  all  contrib¬ 
uting  vibrational  transitions  were  determined  by  a  normalization 
procedure  to  be  described  subsequently. 

(2)  On  the  basis  of  (1),  the  integrated  Intensities  of 
all  contributing  vibration- rotation  lines  were  calculated  system- 
matlcally. 


(3)  As  the  Intensity  and  frequency  of  each  line  was 
generated,  the  line  was  assigned  to  the  frequency  interval  spanning 
it. 


(4)  Once  assigned  to  an  interval,  a  given  line  was  then 
grouped  according  to  the  intensity  decade  into  which  it  fell.  Within 
a  particular  decade  a  running  count  was  made  of  the  number  of  lines 
placed  therein.  Also  retained  was  a  net  sum  of  the  intensities  of 
all  those  lines  placed  in  each  decade.  If  a  line  was  generated  with 
an  intensity  more  than  eight  orders  of  magnitude  below  the  strongest 
line  already  placed  in  the  given  frequency  interval,  it  was  ignored. 
On  the  other  hand,  if  stronger  lines  were  subsequently  generated. 
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the  data  for  weaker  lines  was  dropped  when  they  occurred  In  Intensity 
groups  vAiich  were  more  than  four  groups  below  the  most  intense.  In 
short  either  the  top  five  strongest  groups  were  retained,  or  only 
enough  groups  to  span  eight  orders  of  magnitude.  The  latter  choice 
resulted  in  less  than  five  groups  when  the  lines  in  that  portion  of 
the  spectrum  were  very  weak. 

(5)  The  coefficients  were  then  calculated  for  each  fre¬ 
quency  interval  by  dividing  the  intensity  sums  by  the  number  of  con¬ 
tributing  lines  and  thence  dividing  the  resulting  average  strengths 
by  TXa  Q.  For  H2O  the  value  ■  0.10  cm*l  was  used. 

Once  the  coefficients  ^  and  nj^  were  generated,  the  ef¬ 
fective  transmissions  were  calculated  at  frequencies  corresponding 
to  the  centers  of  the  chosen  frequency  intervals  of  width  5  cm~^. 

(e  ■  0  in  the  notation  of  section  3.2.)  Sequentially  the  coefficients 
1 for  various  combinations  of  the  pressure^  and  amounts  of  absorbing 
gas  were  calculated  (4|(  *  4o)(  u/p).  From  each  of  these  sets  the 
effective  transmissions  were  calculated  as  follows: 

(1)  All  the  transmissions  were  set  equal  to  unity. 

(2)  Starting  with  the  lowest  frequency  Interval,  the  direct 
contributions  to  the  transmissions  were  calculated  at  the  center  of 
the  interval  for  each  intensity  group  contained  therein  using  the 
formulation  described  in  section  3.2  (specifically  Eq.  (32)  and,  if 
applicable  its  limiting  forms,  Eqs.  (33)  and  (35)).  These  were  then 
multiplied  together  and  thence  by  the  factor  1  initially  set.  Hence 
the  resulting  transmission  represented  the  direct  contributions  from 
the  five  intensity  groups  which  represent  all  lines  having  an  intensity 
within  10*^  of  the  strongest  line  in  the  frequency  interval. 

(3)  Subsequent  to  the  calculation  of  the  direct  contribu¬ 
tions,  Che  coefficients  were  used  to  calculate  Che  wing  contributions 
to  the  transmissions  at  the  center  of  the  Interval  of  interest.  The 
wing  contributions  from  lines  at  both  lower  and  higher  frequencies 


It  should  be  mentioned  that  the  term  pressure  is  used  here  in  the 
sense  of  Howard,  Burch  and  Williams,  viz  for  water  vapor  the  effective 
broadening  pressure  is  equal  to  the  sum  of  the  total  pressure  plus 
5.3  times  the  partial  water  vapor  pressure. 
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were  considered.  The  calculation  was  made  from  Eq.  (29)  for  the 
Lorentz  line  shape  and.  In  the  case  of  CO2,  from  Eq.  (42)  for  the 
Benedict  modification.  Thus  as  any  contribution  to  the  transmission 
at  a  given  frequency  was  generated  (be  their  origin  direct  or  wing). 
Its  value  was  multiplied  Into  those  previously  calculated  at  that 
frequency. 


(4)  The  wing  calculations  were  carried  out  at  frequencies 
farther  and  farther  from  the  chosen  Interval  until  their  contribu¬ 
tions  to  the  transmissions  exceeded  0.999. 

Insofar  as  the  normalization  procedure  was  concerned,  the 
detail  yielded  by  the  quasi-random  model  was  not  needed.  This  Is 
because  the  normalization  was  based  on  a  comparison  of  the  calculated 
Integrated  absorptions  with  the  measured  Integrated  absorptions  of 
Howard,  Burch,  and  Williams  for  large  frequency  Intervals  (>  400  cm~^). 
The  correct  normalization  was  obtained  via  a  two  step  process.  A 
large  mesh  size  (20  cm~^)  was  chosen  to  generate  a  separate  coef¬ 
ficient  table  (similar  to  those  Indicated  by  Appendix  A)  for  each  of 
the  seven  transition  groups  In  H2O  and  the  twenty- two  In  CO2  using 
the  relative  Intensities  as  an  Initial  starting  point.  Thence 
Integrated  transmissions  were  calculated  (using  the  simplifications 
just  mentioned)  at  various  effective  pressures  and  paths  corresponding 
to  experimental  results  of  Howard,  Burch  and  Williams  (T  was  always 
chosen  at  300°K) .  These  were  compared  with  the  measured  Integrated 
transmissions  and  the  squares  of  the  deviations  of  ten  to  twenty  such 
sets  were  summed.  The  coefficients ^ multiplied  by  an 
arbitrary  factor  and  the  sum  of  the  squares  of  the  deviations  were 
again  calculated.  Other  multiplicative  factors  were  subsequently 
chosen  and  plotted  as  a  function  of  the  corresponding  sums  of  the 
squares  of  the  deviations.  The  absolute  minimum  was  then  calculated 
graphically  and  a  first  approximation  to  the  normalization  constants 
Cq  (300°K)  deduced  therefrom  for  each  group.  Using  these  first 
approximations,  the  procedure  was  then  repeated  using  a  finer  mesh 
(5  cm~^).  Again  a  least  squares  fit  yielded  the  best  coefficients, 
which  are  listed  In  Table  4.  In  the  event  that  two  groups  overlapped 
a  two-dimensional  generalization  of  the  above  least  squares  method 
was  used. 
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TABLE  4 


WATER  VAPOR  NORMALIZATION  CONSTANTS 


Spectral  Region  C  (300°K) 


6.3/i 

888 

3.2  /i 

16.5 

2.7  M 

1950 

1.87/Lt 

203 

1.38/i 

47.5 

l.lM 

1.87 

.94^ 

1.11 
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SECTION  5 

USE  OF  APPENDICES  AND  INTERPOLATION  PROCEDURES 


5.1  USE  OF  APPENDICES 


The  five  appendices  following  this  section  contain  the 
results  of  the  present  calculation.  Appendix  A  contains  the  co¬ 
efficients  ^  g  (defined  in  Eq.  (28))  and  n  the  number  of  lines  In 
the  Intensity  decade  for  the  top  five  Intensity  groups  of  each 
frequency  Interval  as  a  function  of  frequency.  A  particular  fre¬ 
quency  entry  corresponds  to  the  center  of  an  Interval  extending 
2.5  cm’^  In  both  directions  from  the  center.  Two  sets  of  co¬ 
efficients  are  Included  for  each  temperature  and  correspond  to  two 
meshes  displaced  2.5  cm*^  from  each  other.  For  Identification 
purposes  the  meshes  with  a  center  wave  number  beginning  at  1000  cm~^ 
have  been  labeled  "unshlfted"  \dille  those  starting  at  1002.5  cm*^ 
have  been  labeled  "shifted".  As  before  the  entries  correspond  to 
the  frequencies  at  the  center  of  two  mesh  points  (e.g.  the  entries 
at  1997.5  cm~^  were  calculated  from  the  spectral  lines  lying  between 
1995.0  and  2000.0  cm"^,  etc.).  Note  that  beyond  3400  cm"^  the  co¬ 
efficient  tables  should  be  used  only  with  extreme  caution  unless 
the  values  are  averaged  over  a  wider  frequency  Interval.  Beyond 
3400  cm~^  the  entries  were  calculated  exclusively  on  the  basis  of 
the  symnetrlc  rotator  model  (see  section  2.2).  These  entries  may 
be  used  to  obtain  transmissions  which  are  subsequently  averaged  over 
wider  frequency  Intervals,  itowever,  the  entries  In  themselves  can¬ 
not  be  too  significant  for  Intervals  of  5  cm~^  since  the  positions 
of  the  lines  on  the  symmetric  rotator  model  may  be  in  error  by  more 
than  this  amount.  Nevertheless  entries  beyond  3400  cm*^  have  been 
Included  for  completeness  since  they  formed  the  basis  for  the 
tabulated  average  transmission  values  (Appendices  C,  D,  and  E). 
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Appendix  B  lists,  as  a  function  of  frequency  and  path  thick¬ 
ness  (u.  In  preclpltable  centimeters)  for  three  temperatures  at  seven 
different  pressures,  the  average  transmissions  for  5  cm~^  Intervals 
calculated  on  the  basis  of  the  quasi-random  model.  Note  that  the 
entries  terminate  at  3400  cm~l-,  since  accurate  frequencies  for  the 
stronger  H2O  spectral  lines  were  only  available  up  to  this  frequency 
and  furthermore  that  they  are  listed  every  2.5  cm'^. 


^pendlx  C  lists  transmissions  averaged  over  20  cm'^  at 
Intervals  of  10  cm~^  for  all  the  temperatures  and  pressures  of 
Appendix  B.  However,  the  entries  go  out  to  10000.0  cm~^.  The  averag¬ 
ing  has  minimized  the  uncertainties  In  the  frequencies  of  Individual 
lines  which  Is  Introduced  by  the  symmetric  rotator  description.  In 
effect  the  entries  In  Appendix  C  represent  the  average  transmissions 
that  a  detector  with  an  acceptance  window  of  20  cm~^  and  unit  response 
would  record  If  centered  on  the  tabulated  frequencies. 


Appendices  D  and  E  are  Identical  to  Appendix  C  with  the  excep¬ 
tion  of  the  size  of  Interval  averaged  over.  Appendix  D  lists  entries 
averaged  over  50  cm’^  while  E  lists  entries  averaged  over  100  cm'^. 

In  these  latter  appendices  entries  occur  every  25  cm~^  and  50  cm~^ 
respectively. 


The  experimental  measurements  of  Burch  et  al  for  the  2.7|i 
and  6.3)1  bands  of  water  vapor  are  compared  with  the  theoretical  cal¬ 
culations  for  two  particular  sets  of  values  of  path  length  and  pressure 
In  Figs.  5-1  and  5-2.  The  theoretical  values  plotted  are  averaged 
over  20  cm”l. 


5 . 2  INTERPOLATION  PROCEDURES 

The  extensive  transmission  tables  given  In  the  appendices 
provide  a  suitable  basis  for  interpolation  In  terms  of  pressure, 
temperature  and  path  length.  Interpolation  In  these  parameters  Is 
readily  achieved  via  analytic  methods. 

The  simplest  analytical  procedure  is  that  consisting  of  logarith¬ 
mic  Interpolation  between  non-zero  values  of  the  absorption  A.  Express 
the  absorption  In  the  form 


A  -  1  -  T 


^-f(p,T,u) 


(60) 


vdiere  f  Is  usually  a  complicated  function  of  pressure  (p),  temperature  (T), 


5-2 


5-3 


FIGURE  5-1.  ABSORPTANCE  V  2.7 JX  H^O  BAND 


PERCENT 


FIGURE  5-2.  ABSORPIANCB  OF  b.ifJi  BAND 


5-4 


and  path  length  (u) .  At  points  where  one  or  both  of  the  tabular  values 
of  T  are  identically  equal  to  unity,  linear  interpolation  on  t  pro¬ 
vides  more  than  sufficient  accuracy.  Otherwise  it  is  far  more  accurate 
to  interpolate  f(p,T,u).  This  is  true  because  the  main  part  of  the 
functional  dependence  of  the  absorption  on  p,  T,  and  u  is  usually 
represented  by  the  exponential  factor.  As  only  three  significant 
figures  are  given  in  the  tables,  a  first  order  interpolation  in 
general  provides  sufficient  accuracy.  It  should  be  noted  that  the 
entries  for  varying  pressure  and  path  length  are  given  in  terms  of 
values  that  are  nearly  equidistant  when  plotted  on  a  logarithmic  scale 
(e.g.  1.0,  0.5,  0.2,  0.1,  0.05,  ...). 

From  Eq.  (60)  one  has 


f(P.T,u)  -  -In  A  . 


(61) 


When  the  absorptance  at  constant  p  is  plotted  against  u  on  a  log- log 
graph,  the  resulting  curve  has  only  a  small  curvature  and  displays 
long  linear  sections.  Similar  remarks  apply  to  an  absorptance  curve  at 
constant  u  plotted  against  p.  For  these  reasons,  the  most  accurate 
interpolation  is  obtained  by  taking  equal  logarithmic  Intervals  for 
p  and  u  so  that 

Af  -  AX  +  at  +  Ay  (62) 

where  x  ■  In  p  and  y  «  In  u.  Thus 


f(p  +  Ap,  T+AT,  u+  Au)"  f(p,T,u)  +  A  f  (63) 


and 


T  (p  +  A  p,  T  +  AT,  u  +  Au)  “  1  -  A(p  +  Ap,  T  +  AT,  u+Au)  - 

1  -  exp  -  (f  +  Af)  (64) 
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The  following  examples  illustrate  the  method  described  above. 

a)  Given;  p  ■  1  atmosphere,  T  ■  300°K 

at  u  -  0.1,  T  -  0.583,  A  -  0.417 
at  u  -  0.2,  T  -  0.439,  A  ■  0.561 

Find;  T  at  u  =  0.17  for  p  ■  1  atm.,  T  ■  300°K 


Procedure;  A(O.l) 
A(0.2) 
In(O.l) 
ln(0.17) 
ln(0.2) 


0.417;  f(O.l)  =  .8747 

0.561;  f(0.2)  -  .5780 

-2.3026 

-1.7720 

-1.6094 


Hence ;  Af 


S  f 

S(ln  u) 


A  (In  u) 


0.5780  -  0.8747 
-1.6094  +  2.3026 


(-1.7720  + 


-0.2967 

0.6932 


0.5306 


-0.2271 


Therefore;  f  -  f(O.l)  +  Af  -  0.8747  -  0.2271  -  0.6476 
t(0.17)  -  1  -  -  0.477 


b)  Given;  p  ■  1  atmosphere,  T  ■  300°K 
at  u  -  2.0,  T  -  0.648,  A  -  0.352 

at  u  =  5.0,  T  «  0.425,  A  -  0.575 

p  ■  0.5  atmospheres,  T  ■  300°K 
at  u  -  2.0,  T  -  0.715,  A  -  0.285 

at  u  -  5.0,  T  =  0.545,  A  »  0.455 

Find;  T(p,u)  at  u  =  3.2,  p  =  0.75,  T  »  300°K 

Procedure;  A(l,2.0)  -  0.352;  f  -  1.0441 
A(l,5.0)  -  0.575;  f  -  0.5534 
A(0.5,2.0)  -  0.285;  f  -  1.2553 
A(0.5,5.0)  -  0.455;  f  -  0.7875 
ln(2.0)  -  0.6931;  In(l.O)  -  0.0000 
ln(3.2)  -  1.1631;  ln(0.75)  -  -0.2877 
ln(5.0)  -  1.6094;  ln(0.5)  -  -0.6931 


2.3026) 


5-6 


Hence:  A  f 


A(ln  p) 


•  “>  *  5(1^ 

0.5534  -  1.0441 


hi 

In  I 

(1.1631  -  0.6931) 


1.6094  -  0.6931 

1.2553  -  1.0441  ._o  ,877^ 

^  -0.6931  +  0.0000 

-0.2517  +  0.0877  -  -0.1640 


Therefore:  f(0.75,3.2)  -  f(l,2.0)  +  Af  -  1.0441  -  0.1640  -  0.8801 
T  (0.75,3.2)  -  1  -  -  0.585 


The  data  contained  In  the  Appendices  of  this  report  Is  avail¬ 
able  on  tape.  Inquiries  as  to  the  availability  of  taped  copies  should 
be  sent  to  Gilbert  N.  Plass,  Aeronutronlc  Division  of  Ford  Motor 
Company,  Newport  Beach,  California. 
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APPENDIX  A 

THE  COEFFICIENTS  i  AND  n 
’  o 

FOR  THE  TOP  FIVE  INTENSITY  DECADES 
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COEFFICIENTS 
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A-11 
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COEFFICIENTS  AIV^  FOR  WATER  VAPOR 
T.25CrK  (SHIFTED  MESH) 
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A-14 


TED 
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APPENDIX  B 

WATER  VAPOR  TRANSMISSIONS  FOR  FREQUENCIES 
LESS  THAN  3400  cm'^  AT  300°K,  250°K,  AND  200°K 
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WATER  VAPOR  1.00  ATMOSPHERE  PRESSURE 
T=300"K 
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WATER  VAPOR  AT  0.50  ATMOSPHERE  PRESSURE 
T=300-K 
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WATER  VAPOR  AT  0.50  ATMOSPHERE  PRESSURE 
T=300-  K 


WATER  VAPOR  AT  0.20  ATMOSPHERE  PRESSURE 
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WATER  VAPOR  AT  0.20  ATMOSPHERE  PRESSURE 
T=300-  K 
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WATER  VAPOR  AT  0.02  ATMOSPHERE  PRESSURE 
T=300‘K 
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WATER  VAPOR  0.01  ATMOSPHERE  PRESSURE 
T=30CrK 
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WATER  VAPOR  0.01  ATMOSPHERE  PRESSURE 
T=300’  K 
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WATER  VAPOR  AT  1.00  ATMOSPHERE  PRESSURE 
T=250”  K 
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WATER  VAPOR  AT  1.00  ATMOSPHERE  PRESSURE 
T  =  250°  K 
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WATER  VAPOR  AT  0.02  ATMOSPHERE  PRESSURE 
T-20Gf  K 
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WATER  VAPOR  AT  0.C2  ATJ 
T-20a’  K 


^  !!’!! !!'!!  Will !-!!!  I!!-  ^ 

i  iiiii  iiiii  iiiii  !!i!!  inn  iiiii  iii!!  iliii  ii!!!  iHii  iiii!  iii!!  wm  iiii!  Hiii  Iw  iii!!  iii!!  iiii!  ww.  !!!!*  ijiii  iilii  iiiii 

i  Villi  niis  mil  silil  riili  iflii  llm  Ifiii  iiiii  HiH  iiiii  iHii  iiiii  iiiii  iiiH  iiiii  iHI;  iiiii  niil,  iiiii  iiiii  iiiii  liiii  {Mli  iiii; 

^oeea  «e««^  eaeao  aaaa^  Jaaaa  aa^aa  aaaaa  a^j4a  aaaaa  .aava  •  .J  — _Jaaa  a----  -aaaa  HHm  --JJii  JlmiSi  mmHo  aJJ>-  aaaa>  J^aJa 


WATER  VAPOR  AT  0.02  ATMOSPHERE  PRESSURE 
T-200r  K 
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WATER  VAPOR  AT  O.OI  ATMOSPHERE  PRESSURE 
T-200-  K 
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WATER  VAPOR  AT  0.01  ATMOSPHERE  PRESSURE 
T-20Cr  K 
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APPENDIX  C 

WATER  VAPOR  TRANSMISSIONS  AVERAGED  OVER  20  cm"’' 
WITH  ENTRIES  EVERY  10  cm"’'  FOR  FREQUENCIES  FROM 
1000  cm"’'  TO  10,000  cm"’'  AT  300°K,  250°K,  AND  200°K 
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WATER  VAPOR  AT  0.50  ATMOSPHERE  PRESSURE 
T=30(r  K 


WATER  VAPOR  AT  0.20  ATMOSPHERE  PRESSURE 
T=300-  K 
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WATER  VAPOR  0.10  ATMOSPHERE  PRESSURE 

T=300’K  AVERAGED  OVER' 20cm 
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WATER  VAPOR  AT  0.10  ATMOSPHERE  PRESSURE 

T=30(rK  AVERAGED  OVER'  20cm 
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WATER  VAPOR  AT  0.05  ATMOSPHERE  PRESSURE 
T=30CrK 
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AVERAGED  OVER  20cm' 


WATER  VAPOR  AT  0.02  ATMOSPHERE  PRESSURE 
T=300'  K 
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WATER  VAPOR  AT  0.01  ATMOSPHERE  PRESSURE 
T=300r  K 
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C-23 


WATER  VAPOR  AT  1.00  ATMOSPHERE  PRESSURE 
T-25a’  K 
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WATER  VAPOR  AT  I.OO  ATMOSPHERE  PRESSURE 
T  =  25Cr  K 
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WATER  VAPOR  AT  1.00  ATMOSPHERE  PRESSURE 

T-250’  K  AVERAGED  OVER  20cin 
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WATER  VAPOR  AT  0.10  ATMOSPHERE  PRESSURE 
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WATER  VAPOR  AT  0.02  ATMOSPHERE  PRESSURE 
T  =  200°  K 
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WATER  VAPOR  AT  0.02  ATMOSPHERE  PRESSURE 
T-20a  K 
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WATER  VAPOR  AT  0.02  ATMOSPHERE  PRESSURE 
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APPENDIX  D 


WATER  VAPOR  TRANSMISSIONS  AVERAGED  OVER  50  cm 
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WITH  ENTRIES  EVERY  25  cm  FOR  FREQUENCIES  FROM 
1000  cm‘^  TO  10,000  cm’^  AT  300°K,  250°K,  AND  200°K 
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WATER  VAPOR  AT  100  ATMOSPHERE  PRESSURE 
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WATER  VAPOR  AT  0.20  ATMOSPHERE  PRESSURE 
T=300"  K 
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WATER  VAPOR  AT  0.02  ATMOSPHERE  PRESSURE 
T=300’  K 


WATER  VAPOR  AT  0.01  ATMOSPHERE  PRESSURE 
T=30(rK 
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WATER  VAPOR  0.50  ATMOSPHERE  PRESSURE 
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WATER  VAPOR  AT  0.20  ATMOSPHERE  PRESSURE 
T=2S0P  K 
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WATER  VAPOR  AT  0.10  ATMOSPHERE  PRESSURE 
T=250“  K 
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WATER  VAPOR  AT  0.10  ATMOSPHERE  PRESSURE 

T-200^K  cr,  ■ 

AVERAGED  OVER  50cm 
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WATER  VAPOR  AT  0.05  ATMOSPHERE  PRESSURE 
T-200f  K 
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WATER  VAPOR  AT  O.Ot  ATMOSPHERE  PRESSURE 

T=200°  K  AVERAGED  OVER  50cm  ' 
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APPENDIX  E 

WATER  VAPOR  TRANSMISSIONS  AVERAGED  OVER  100  cm'^ 
WITH  ENTRIES  EVERY  50  cm‘^  FOR  FREQUENCIES  FROM 
1000  cm"^  TO  10,000  cm'^  AT  300°K,  250°K,  AND  200°K 
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WATER  VAPOR  0.10  ATMOSPHERE  PRESSURE 

T^SOOfK  AVERAGED  OVERlOOcm  ’ 
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WATER  VAPOR  AT  0.05  ATMOSPHERE  PRESSURE 
T=30Cr  K 


WATER  VAPOR  0.01  ATMOSPHERE  PRESSURE 

T=300"K  AVERAGED  OVER  100cm 
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water  vapor  at  >00  atmosphere  pressure 

T-25Cr  K 
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WATER  VAPOR  AT  0.50  ATMOSPHERE  PRESSURE 
T-250’  K 
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WATER  VAPOR  AT  0.20  ATMOSPHERE  PRESSURE 
T=250°K 
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WATER  VAPOR  AT  0.10  ATMOSPHERE  PRESSURE 
T=250°  K 
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WATER  VAPOR  AT  0.02  ATMOSPHERE  PRESSURE 
T=250°  K 
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water  vapor  at  O-'O  atmosphere  pressure 
T=20af  K 
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Absorption 
Water  Vapor 
Infrared  Spectroscopy 


Final  Report,  Volume  11,  20  September  1962 

249  p.,  incl.  tables  and  ILlus.  Unclassified  Report 


The  infrared  absorption  of  water  vapor  has  been  calculated 
over  a  wide  range  of  path  length,  pressure,  and  temperature 
from  1000  to  10,000  cm"^.  The  results  are  presented  in 
extensive  tables  which  give  the  transmission  over  intervals 
of  2.S  cm*^  at  wave  numbers  below  3400  cm~^.  the  average 
transmission  over  intervals  of  20,  SO,  and  100  cm~^  is  given 
over  the  complete  frequency  range.  The  calculated  data 
representing  the  intensity,  position,  and  number  of  spectral 
lines  in  a  given  frequency  Interval  is  tabulated  in  a  form 
which  is  useful  for  transmission  calculations.  All  lines 
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from  all  isotopic  species  having  an  intensity  greater  than 
10*^  of  Che  strongest  line  in  an  absorption  region  were 
included  in  the  calculacion.  The  quasi^random  model  Cor  the 
representation  of  band  spectra  is  discussed  in  detail.  Among 
Its  features  are  the  possibility  of  representing  the  absorp¬ 
tion  by  a  series  of  lines  which  are  arranged  neither  at 
random  nor  in  a  regular  way. 
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